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Abstract 

We explain Barkhausen noise in magnetic systems in terms of avalanches 

near a plain old critical point in the hysteretic zero-temperature random-field 

Ising model. The avalanche size distribution has a universal scaling function, 

making non-trivial predictions of the shape of the distribution up to 50% 

above the critical point, where two decades of scaling are still observed. We 

simulate systems with up to 1000 3 domains, extract critical exponents in 2, 3, 

4, and 5 dimensions, compare with our 2d and 6 — e predictions, and compare 

to a variety of experimental Barkhausen measurements. 
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When materials are pushed, they can yield in different ways. Some crackle: they trans- 
form through a series of pulses or avalanches. In many systems, the behavior of these 
avalanches is unaffected by thermal fluctuations: one domain triggers, pushing some of its 
neighbors to trigger, in a deterministic process dependent on the static, quenched disorder in 
the material (and on the stress history). The statistics of Barkhausen noise (the avalanches 
seen in magnetic materials as the external magnetic field is ramped up and down) has 
been extensively studied experimentally |I|,3,^^^|6|,|7|,|^J^ JlO|JTl|| . We suggest that the zero- 



temperature random-field Ising model |12||L3| provides a universal, quantitative explanation 
for many of these experiments. 

A typical experiment will collect a histogram of pulse sizes, times, or energies. The 
distribution will follow a power law, which cuts off after two to several decades — much 
broader than any observed morphological feature in the materials. An explanation for the 
experiment must involve collective motion of many domains; it must provide an explanation 
for the observed power-law scaling regions, and it must provide an explanation for the cutoff. 

Figure 1 shows the distribution D int (S,R) of avalanche sizes for our model in 3d (dis- 
cussed already in reference [|12||), at several values of the microscopic disorder R. The model 
is a collection of domains Sj = ±1 coupled to an external field H, a local random field fa 
chosen from a distribution p(f) = exp(— f 2 /2R 2 ) /\/2iiR of standard deviation R, and to its 
nearest neighbors Sj with an energy of strength J = 1. The domain s» flips over when the 
net local field Fj = fi + H + JJ2 nn s j seen a t site i changes sign. Due to the nearest-neighbor 
interaction, a flipping spin often causes one or more neighbors to flip also, thereby spawning 
a whole avalanche of spin flips. Figure 1 shows the avalanches found by integrating as the 
external field H(t) is raised adiabatically from — oo to oo (the field is thus constant during 
the individual avalanches). 

Notice three things about figure 1. (1) The distributions follows a power law, which 
cuts off after two to several decades. (2) The cutoff appears to diverge at a critical value 
of the disorder R c , which we estimate in three dimensions to be 2.16J. (3) The critical 
region is large! While the true power-law distribution is only obtained at R c = 2.16, we get 



avalanches with more than a hundred domains all the way up at R = 4. This suggests that 
experiments can see decades of scaling without working hard to find the critical disorder. 
Several decades of scaling without tuning a parameter need not be self-organized criticality: 
it can be vague proximity to a plain old critical point. 

Notice four more things about figure 1. (4) The straight line lying askew below the 
numerical data is our prediction for the asymptotic power law Di nt (S, R c ) ~ s -( r + cr / 35 ). The 
obvious experimental method of taking the slope on the log-log plot gives the wrong answer 
until many, many decades of scaling are obtained! (5) The inset shows the collapse of the 
data onto a scaling function 

V mt {s°{R - Rc)/R) = Jim s T+ ^ s D mt (S, R) (1) 

it — *Rc 

which is a universal prediction of our model: real experiments rescaled in the same way 
should look the same (apart from overall vertical and horizontal scales). This scaling function 
is quite unusual: it grows by a factor of over ten before cutting off. This bump is the reason 
that the experiments take so long to converge to the asymptotic power law. To make a 
definite prediction, we have phenomenologically fit our curve 

V int (x) = (0.021 + 0.002x + 0.531x 2 - 0.266x 3 (2) 
+ 0.261x 4 ) exp(-0.789x 1/CT ) 

where we guess the error in the curve to be less than 10% within the range 0.2 < x < 1.2. 
(6) In the main figure, the scaling form passes through our data quite well, even far from 
R c . The scaling theory is predictive for curves with only two or three decades of scaling. 
The critical region starts when the correlation length (and hence the avalanche size cutoff) 
becomes large — not only when the pure power law takes over. Using equations ([I]) and (2) 
and the values a = 0.24 ± 0.02 and r + a/35 = 2.03 ± 0.03, an experimentalist should be able 
to fit any single histogram of avalanches, shifting only the overall vertical and horizontal 
scale factors. (7) Widom scaling (equation (]!])) forms a powerful tool: only by varying R 
were we able to extract the correct critical behavior. We suggest that experimentalists try 
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varying some parameter of their system (annealing time or temperature, grain size, impurity 
concentration, ...) and observe the resulting cutoff dependence. Any family of curves thus 
generated should, near the critical point, be fit with three parameters (including R c in (1) 
and (2)). 

A comparison of our predicted exponents with power laws extracted from a number of 
experiments on magnetic Barkhausen noise in bulk three-dimensional systems is shown in 
Table 1. One of the experiments |§ varied the annealing time, and saw a shift in the cutoff, 
but did not extract a critical annealing time or do scaling collapses. The range of values 
for r + a/35 observed in the experiments is compatible with the range of log-log slopes we 
observe due to the unusual scaling form for the integrated avalanche size distribution D int 
discussed above. 

The largest set of experiments measure the avalanche size distribution in a narrow range 
of fields (i.e., without averaging over the entire loop): their power laws fits are a measure of 
the critical exponent r. Integrating over the hysteresis loop changes the power law, because 
only near a critical value H c of the external field do large avalanches occur. The cutoff in the 
avalanche size at R c goes as | H — H^ 1 ^ 135 : a as discussed above gives the cutoff dependence 
with R — R c , and (3 and 1/5, as usual |T2|], give the singularities of the magnetization with 



R — R c and H — H c respectively (the exponent for any quantity varying H at R c is given by 
multiplying 1/(35 by the exponent for the singularity varying R at H c ). The scaling function 
for the non-integrated avalanche size distribution ||12|| , we find, does not have a bump. The 



experimental measurements for r are close to our numerical estimate. 

The other experiments (pulse durations, power spectra, and pulse energies) introduce a 
new exponent combination avz. The correlation length exponent v gives the divergence of 
the characteristic spatial extent of avalanches with R — R c , and vz gives the divergence of 
the avalanche durations with R — R c . The critical exponent z occasionally can depend on the 
details of the spin dynamics [15[]; it is not even clear whether our simulation must have the 



same value of z as our e-expansion. Nonetheless, the agreement between our predictions and 
the measured values are about as good as the agreement between the different measurements. 



We have also investigated the application of our model to other systems. Many experi- 
ments are done in effectively two-dimensional systems (magnetic hysteresis ||16|| , avalanches 
as the field is swept in superconductors ||17|| , and avalanches as helium is injected into Nu- 
clepore |18[); our 2d explorations are still rather preliminary. Our model does not fit the 
avalanche size distributions measured in 3d martensitic transitions as the temperature is 
ramped |19|]: their measurement of the pulse duration exponent (r + a(35 — l)/avz + 1 ~ 1.6 
is significantly different from our prediction of 2.81 ± 0.11. We expect that the long-range 
anisotropic elastic fields in martensites likely change the universality class; similarly the 
long-range antiferromagnetic demagnetization fields could affect experiments in certain ge- 
ometries (see ]l(|). Full explanations about the various exponent combinations measured 
in different systems |20| and detailed discussions of experiments and systems [21| are forth- 



coming. 

Figure 2 shows the results for five of our exponents (r + o (35, r, 1/z/, av, and avz), 
measured in 2, 3, 4, and 5 dimensions. (From these one can get (3 and 5 separately using the 
exponent equality |^U| r + a {35 = 2 + a (3.) We measure the exponents in the (unphysical) 



dimensions four and five in order to test our renormalization-group predictions [ ]I3"| , PU] for 
the behavior near six dimensions. First, notice that the numerical values converge nicely 
to the mean-field predictions as the dimension approaches six, and that the predictions of 
our 6 — e expansion do remarkably well. (The primary role of the renormalization-group 
treatment, of course, is to explain why scaling and universality might be expected in these 
systems.) The predictions for 1/z/ are to fifth order in e: by mapping our model to all orders 
onto the regular Ising model in two lower dimensions |22|] , we have been able to use [23 



the series known to order e 5 for v. The other exponents shown have no equivalent in the 
equilibrium model: we have developed [^(J a new method for calculating these avalanche 
exponents using two replicas of the system. The dashed lines show a Borel resummation 
f23| of the series for 1/z/, and the predictions to first order for the other variables. 

Second, notice the exponents in two dimensions. We here conjecture that the 2d expo- 
nents will be t + a (35 = 2, r = 3/2, 1/z/ = 0, and av = 1/2. It is likely that two is the 



"lower critical dimension" for our system, below which all avalanches will be finite except at 
zero disorder. At the lower critical dimension, the critical exponents are often ratios of small 
integers, and it is often possible to derive exact solutions. For us, using the fact that there 
can be at most one system-spanning avalanche in two dimensions, one can derive a special 
exponent relation 1/av — d — (3/v = 2 — (3 jv: this "hyperscaling" relation is false in mean- 
field theory, and definitely ruled out numerically in four and five dimensions, and probably in 
three p4[,|l4f . Folklore in the field |25[ give us two other likely 2d relations: one each for the 
exponent giving the decay in space of the cluster correlation function fj = fl/v+4—d =? A—d 
and of the avalanche correlation function 1] = 2 + /3/v — (35 jv =? 1. These relations hold in 
the lower critical dimension for the Ising model, the Heisenberg model, and the equilibrium, 
thermal random-field Ising model. The first of these relations is equivalent to the statement 
that the avalanches are compact (l/au = d = 2). 

We must mention that our firm conjectures about the exponents in two dimensions must 
be contrasted with our lack of knowledge about the proper scaling forms. At the lower 
critical dimension, the correlation length typically diverges exponentially as one approaches 
the critical point (v — > oo). Some combinations of critical exponents stay finite (hence 
au = 1/2), but those which diverge and those which go to zero usually must be replaced 
by exponents and logs, respectively. We have used three different RG-scaling ansatze to 
model the data in two dimensions. (1) We used the traditional scaling form £ ~ (R — R c )~ u , 
deriving v = 5.3 ± 1.4 and R c = 0.54 ± 0.04. These collapses worked as well as any, but the 
standard form has more free parameters to fit with. Also, the large value for v (and larger 
values still for I /a = 10 ± 2) makes one suspicious. (2) We used a scaling form suggested 
by Bray and Moore [[26[ in the context of the equilibrium thermal random-field Ising model 



where R c = 0: if they assume that R is a marginal direction, then by symmetry the flows 
must start with i? 3 , leading to £ ~ exp(A/(R— R c ) 2 ) = exp(A/ R 2 ). This form had the fewest 
free parameters, and most of the collapses were about as good as the others (except notably 
for the finite-size scaling of the moments of the avalanche size distribution, which did not 
collapse well once spanning avalanches became common). (3) We developed another possible 
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scaling form, based on a finite R c and R marginal, which generically has a quadratic flow 
under coarse-graining: here £ ~ exp(A/ (R — R c ). Here again the moments did not scale well; 
we find R c = 0.54 ± 0.04, quite compatible with the traditional scaling collapse. This is not 
a surprise: it is always hard to distinguish large power laws from exponentials. Amazingly 
enough, the exponents plotted in figure 2 were largely independent of which scaling form we 
used! The error bars shown span all three ansatze, and are compatible with our conjectures 
above. 

We are not the only ones to model avalanche behavior in disordered magnets. There has 
been much work on depinning transitions and the motion of individual interfaces p7| , |28| ; 
our system, with many interacting interfaces, perversely seems much simpler to analyze. 



Many have studied related models with random bonds p9|j30|| and random anisotropies; 
random fields are actually rather rare in experimental systems. We now believe on symmetry 
grounds that all these systems are in the same universality class (as argued numerically 



and previously shown for depinning pq|). The external field H c at the critical point breaks 
the rotational and up-down symmetries of these models (and of the experiments!), and the 
spins which flip far from the critical point (roughly M(H C )) act as random fields. On the 
other hand, we ignore long-range forces (discussed above) and long-range correlations in the 
disorder (e.g., dislocation lines and grain boundaries): these likely will lead to closely related 
but distinct universality classes. 

We acknowledge the support of DOE Grant #DE-FG02-88-ER45364 and NSF 
Grant #DMR-91 18065. We would like to thank Bruce W. Roberts, Sivan Kartha, 
Eugene Kolomeisky, James A. Krumhansl, Mark Newman, Jean Souletie and Uwe 
Tauber for helpful conversations. This work was conducted on the SP1 and the 
SP2 at the Cornell National Supercomputer Facility, funded in part by NSF, by 
NY State, and by IBM. Further pedagogical information using Mosaic is available at 
http : / / www. lassp . Cornell. edu / set hna / hysteresis / hysteresis . html . 
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TABLES 



exponents 


simulation 


experiments 




in 3 dim. 


in 3 dimensions 



T + (7/3(5 

(r - l)/avz + 1 

(r + a(3S- l)/auz + 1 
(3 — t)/ctisz 
(3- {r + aPS))/avz 
(r - l)/(2 - auz) + 1 



1.6 ±0.06 



2.03 ±0.03 
2.05 ±0.12 

2.81 ±0.11 
2.46 ±0.17 
1.70 ±0.10 
1.42 ±0.04 



1.74,1.78,1.88 @: 

1.5 ±0.5 1; 
1.33 ||; 
1.5- 1.7 § 
1.73-2.1 § 
1.64,2.1,1.82 g; 
1.7-2 § 

2.28 § 
around 2 |J| 

1.6 §§; 1.8 
1.44,1.58,1.60 § 



TABLE I. Critical exponents obtained from numerical simulations [12,14] and ex- 
periments on Barkhausen noise in different magnetic materials (Fe, alumel, metglass [Q], NiS [||, 
SiFe |,§|7|, 81%NiFe §, AlSiFe g, and FeNiCo The sample shapes were mostly wires. The 

quoted exponents were experimentally obtained from the pulse-area distribution in a small bin of 
the magnetic field H (exponent r), the pulse-area distribution integrated over the entire hysteresis 
loop (r + a(35), the distribution of pulse durations in a small bin of H ((r — l)/avz + 1), the dis- 
tribution of pulse durations integrated over the loop ((r + af35 — l)/avz + 1), the power spectrum 
of the pulses in a small bin of H ((3 — t)/(jvz), the power spectrum of the pulses integrated over 
the hysteresis loop ((3 — (t + a (35)) /auz), and the distribution of pulse energies in a small bin of 
H ((r — l)/(2 — crvz) + 1). Notice that these experiments are mostly done in geometries which 
minimize the effects of demagnetization fields (deadly to our model). 
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FIGURES 

FIG. 1. Avalanche size distribution curves in 3 dimensions integrated over the external field. 
From left to right, the first three curves are for system size 320 3 and disorders 4.0, 3.2, and 2.6. 
They are averages over different initial random field configurations. The last curve is a 1000 3 run 
at R = 2.25, where R c = 2.16 and r = (R — R c )/R- The inset shows the scaling collapse of curves 
in 3d. The disorders range between 2.25 and 3.2; the top curves at R = 3.2 show noticeable 10% 
corrections to scaling. In the main figure, the avalanche size distribution curves obtained from the 
fit to this data (thin lines) are plotted alongside the raw data (thick lines). Notice that the scaling 
theory is predictive up to R = 3.2, 50% above R c . The long-dashed straight line is the expected 
asymptotic power-law behavior s~ 2 ' 03 . Notice that it does not agree with the measured slope of 
the raw data. 

FIG. 2. Numerical values (filled symbols) of the exponents r+cr/35, r, ovz, and av (circles, 
diamond, triangles up, squares, and triangle left) in 2, 3, 4, and 5 dimensions. The empty symbols 
are values for these exponents in mean field (dimension 6). Note that the value of r in 2d was not 
measured. The empty diamond represents the conjectured value (see text). We have simulated 
sizes up to 7000 2 , 1000 3 , 80 4 , and 50 5 , where for 320 3 for example, more than 700 different random 
field configurations were measured. The long-dashed lines are the e expansions to first order for 



the exponents r + a (35, r, avz, and av. The short-dashed line is the Borel sum 23] for 1/v to 
fifth order in e. The error bars denote systematic errors in finding the exponents from collapses of 
curves at different values of disorder R. Statistical errors are smaller. 
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